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COMPARISON OF METHODS OF FORECASTING NIGHT 
MINIMUM TEMPERATURES ON CONCRETE ROAD SURFACES 


By G. E. PARREY, W. G. RITCHIE and S. E. VIRGO, O.B.E. 


Summary. Three methods of forecasting night minimum temperatures at the surfaces of 
concrete roads were examined and possible improvements considered. A regression-equation 
method which gives the forecast concrete minimum temperature as a function of the air temp- 
erature and dew-point at 12 Gmr the previous day was improved by making an adjustment for 
cloud amount during the night. Satisfactory improvements were not obtained for the other two 
methods considered which derive the forecast concrete minimum temperature from a combina- 
tion of a forecast of the air minimum temperature with a forecast of the difference between 
the air minimum temperature and the concrete minimum temperature. The methods were 
applied over periods independent of the original data, and successful results were obtained at 
two independent road surfaces. The mean errors and standard deviations of the errors are 
given. The methods depending on forecast of air temperature give better results than the 
direct regression method unless there are large errors in forecasting air temperature. 


Introduction. First reports on two experiments concerned with night 
minimum temperatures at the surfaces of concrete roads have already 
appeared in this magazine; Parrey’ has reported on a road at Watnall and 
Ritchie® has reported on one at Wyton. Both experiments are continuing and, 
as time goes on, more data are accumulating. Since the practical application 
is to forecasting temperatures of o°C or below on roads or runways, this 
paper will be concerned only with the months October to April inclusive. 


In both experiments ordinary minimum thermometers are used with their 
bulbs resting on the surfaces of the roads. Other reports on minimum temp- 
eratures at or near concrete surfaces have also been published, for example 
by Johnson and Davies* and by Hay.* In these experiments temperatures 
were measured by other means, but as it is not yet known whether results 
obtained by the various methods are strictly comparable, this paper will 
confine itself to a discussion of the Wyton and Watnall results. 


There are two ways of attacking the problem of forecasting night minimum 
temperatures over concrete. One is to forecast the concrete minimum 
temperature directly by means of a regression equation from observations 
made in a thermometer screen. The other is to forecast first the air minimum 
temperature by one of the recognized methods and then as a second step to 
derive the concrete minimum temperature from it. 


349 
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The direct regression method. For the first method, the direct method, 
data from Watnall for the period October 1967-April 1968 were examined 
and occasions when a front passed between 12 Gmt and 06 GMT next morning 
were rejected. From the remaining 147 cases the following correlation 
coefficients were obtained : 


between Mz and 7,, 0°75 
between Mz and D,, 0-81, 


where Mz is the minimum temperature at the surface of the concrete road, 
Ty, is the dry-bulb temperature at 12 GMT 
and D,, is the dew-point temperature at 12 GMT. 


These correlation coefficients suggested that a useful regression equation 
could be set up and the following equation was obtained by the method of 
least squares? : 

Mr = 0°59T 2 + 0°69D,2 - 46, 
where temperatures are in degrees Celsius. 


Values derived from this equation were tabulated for Watnall in Table I 
for ready use. Differences between forecasts obtained from this equation and 
observed minimum temperatures were calculated and an attempt was made 
to censtruct a correction table to make some allowance for the effects of cloud 
amount and wind speed. There was no apparent relation with wind speed 
and the residuals were very widely scattered in relation to cloud amount. 
The best that could be done was to add the small correction table to the 
main Table I. (Presumably the reason why no relation with wind speed was 
found is that wind direction is important; if the wind blows along the concrete 
road, the fetch is over concrete; but if it blows across the road, the fetch is 
over grass.) 


TABLE I-—TABLE FOR OBTAINING NIGHT MINIMUM TEMPERATURE AT WATNALL 
AT THE SURFACE OF CONCRETE FROM AIR TEMPERATURE AND DEW-POINT 
TEMPERATURE AT I2 GMT ON THE PREVIOUS DAY 
Air temp- Dew-point temperature D;(°C) 
erature 7) -7 -6 =-5§ -4 -3 -2 -!1 0 1 2 
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Correction (8) to be added for mean cloud amount during the night: 
Mean cloud amount (oktas) 0 1 3 4 5 6 7 8 
8 (degC) -12 -07 00 00 403 +06 +09 +1-1 +13 


The next step was to test this equation on samples of data obtained from 
two separate winters at Watnall and three separate winters at Wyton. In 
every case the cloud correction produced a marginal improvement in the 
standard deviation of the errors, but moved the mean a little further from 
zero; this is not surprising since the basic equation was derived by the method 
of least squares. 
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The winter of 1969-70 was one of those for which data from Watnall were 
used. October 1969 contained many warm days in a range not represented 
in the data used for deriving the equation, so a recalculation of the data for 
this winter was carried out by rejecting occasions when D,, was greater than 
10°C. This was deemed permissible as forecasting temperatures of o°C or 
below on roads is unlikely to be a problem for a forecaster in these circum- 
stances, and moreover the regression may not be strictly linear at these extreme 
values. Table II shows the improvement which resulted. 


TABLE II—EFFECT OF REJECTING OCCASIONS WHEN D,, WAS GREATER THAN 10°C 
IN THE APPLICATION OF THE DIRECT REGRESSION METHOD TO WATNALL DATA 
FOR THE PERIOD OCTOBER 1969—APRIL 1970 


Number of Mean Standard 
observations error deviation 
degrees Celsius 
Without cloud correction : All cases 110 +060 3°54 
Rejecting 
Dy,>10°C — 096 — 0°04 2°83 
With cloud correction : All cases +103 3°43 
Rejecting 
Dy,>10°C —_ 96 +0°29 2°66 


Analysis of variance showed that the three samples of data consisting of 
errors in forecasts of minimum temperatures at Wyton all came from the 
same population. Data for all three winters were then lumped together to 
obtain the best estimates of the mean and standard deviation of this population. 
Analysis of variance also showed that the data for two winters at Watnall 
were themselves samples of one population and the mean and stande:d 
deviation were calculated for this population. In short, there was no year-to- 
year variation at either place. Student’s t-test was now applied to these two 
resultant populations and it was found that these two populations were them- 
selves samples of a common population. The whole process was then repeated 
with allowance for cloud cover and this process again yielded resultant 
populations for Wyton and Watnall which were themselves samples of a 
common population (though not the same population as the first which con- 
tained no correction for cloud cover). The 5 per cent level of significance was 
used as the criterion and, as stated above, occasions when D,, was greater 
than 10°C were rejected. 


It is a noteworthy result that the errors in forecasts of minimum temperatures 
on two roads over 60 miles apart are samples of the same two populations 
(with and without cloud cover respectively); but as later work* showed that 
there are differences between the results for other places in the same part of 
England by this method, this particular result may be no more than a 
coincidence. Throughout this paper the error is reckoned as the forecast 
value minus the observed value, and the best estimates of the mean error and 
standard deviation for each of the two populations are as follows : 


Without cloud With cloud 
correction correction 
Mean error - aes degC +0-38 degC 
Standard deviation 2°50 degC 2°42 degC 
Total number of cases 554. 
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The reduction in standard deviation is the maximum which can be expected 
from applying the cloud correction because cloud amounts were assessed 
from the Daily register after the event and are therefore not subject to any 
forecasting errors. 


Because Parrey and Ritchie had shown that the depression of the minimum 
temperatures at the surface of a concrete road below that of the air in the 
screen is closely related to the hours of darkness, an attempt was made to 
introduce a further modification into the direct regression equation to take 
the duration of darkness into account. It was unsuccessful. Perhaps an 
equation might have been set up to derive Mg in terms of three quantities 
instead of two (duration of darkness perhaps being the third) but an equation 
of this kind cannot be tabulated for ready use on the forecast bench, and the 
figures given above therefore probably represent the best that can be achieved 
by direct regression. 


Indirect methods. (i) Parrey’s method. Two methods have been proposed 
for deriving a forecast of the minimum temperature over a concrete road 
(Mr) from a forecast of the minimum air temperature (M4) in the 
thermometer screen. From data for the winter 1967-68 Parrey' derived the 
linear regression equation : 


Ma—- Mp = 0°28t - 2-9, 


where ¢ is the time between sunset and sunrise in hours and temperatures 
are in degrees Celsius. With the aid of the Nautical Almanac this equation 
can be tabulated for ready use for any given latitude. Table III shows the 
tabulation for 53°N. This method was not tested on the road at Watnall 
but it was tested on data for three successive winters at Wyton. As this was 
done in retrospect, actual values of M, were used. Analysis of variance 
showed that data for all three winters could be regarded as samples of the 
same population and the following results were obtained for all occasions : 


number of occasions 618 
mean error -— 0-04 degC 
standard deviation c, 1:02 degC 
The errors were normally distributed. 


To obtain the total error it is necessary to take into account the errors in 
forecasting air minimum temperatures. Steele, Stroud and Virgo’ have given 
figures for clear and cloudy nights separately for the period October 1967- 
March 1968. When compounded the mean error is — 1-10 degC and the 
standard deviation o, = 2:79 degC. By adding the errors and adding the 
variances we may obtain an estimate of the errors involved in the whole 
process. The figures are : 


mean error — 1°14 degC 
standard deviation 4/(o,? + o,") = 2-97 degC. 


(ii) Ritchie’s method. Instead of a regression equation Ritchie? expressed 
the quantity M,4—- Mp for each day of the year October 1967—September 
1968 as a simple sine curve of the form 


Ma —- Mr = 0°48 + 1-22 sin (6 + 9), 
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TABLE III—TABULATION OF (M,4—- Mr) FOR PARREY’S METHOD 
Duration between sunset Ma- Mr Dates 
and sunrise 

from degC 

h min. 

g 06 

9 27 

9 49 

10 10 


32 


16 


38 
oo 
21 


41 
25 


to 
h min. 
g 26 


26-30 Apr. 
9 48 


20-25 Apr. 
10 09 14-19 Apr. 
10 31 9-13 Apr. 
10 $3 4-8 Apr. 
11 15 29 Mar.—3 Apr. 
II 37 24-28 Mar. 
11 59 19-23 Mar. 
12 20 14-18 Mar. 
40 9-13 Mar. 
03 4-8 Mar. 
24 27 Feb.-3 Mar. 
45 21-26 Feb. 
06 16-20 Feb. 


1-5 Oct. 
6-10 Oct. 
11-16 Oct. 
17-21 Oct. 
22-26 Oct. 


+o°8 
+0o'9 
+10 
+11 


27 
49 
10 
32 
5 
I 


+1°2 
+13 
+14 
+15 
+16 


10-15 Feb. 


4-9 Feb. 
29 Jan.—3 Feb. 
22-28 Jan. 
15-21 Jan. 
8-14 Jan. 


27-31 Oct. 
1-6 Nov. 
7-12 Nov. 
13-19 Nov. 
20-27 Nov. 
28 Nov.-8 Dec. 


4 
16 ae +1°7 9 Dec.—7 Jan. 
The durations on the left-hand side of the table are derived from the equation 
Ma - Mr = 0°28t— 2-9 where ¢ is duration in hours. The corresponding dates on 
the right-hand side are appropriate to latitude 53°N. The quantity in the middle 
column must be subtracted from Mz, to obtain Mr. 


where 6 is the day as an angle (365 days = 360°) and 9 is a phase angle. 
When a harmonic analysis was performed on the data for another year, 
October 1968-September 1969, it was found that the numerical constants 
were the same as before to two places of decimals and the only difference 
was in the phase angle (due presumably to a late spring or early autumn or 
some such difference between the years). This curve, with the value of the 
phase angle obtained for the year October 1967—September 1968, was taken 
as the basis for further study. The actual points for the year October 1967- 
September 1968 are scattered about this curve and an attempt was made to 
reduce the scatter by making allowances for wind speed and cloud amount. 
This was unsuccessful, so the best forecast of the depression of Mp below M4 
by this method remains that derived from the curve itself. Values of this 
quantity were tabulated for 10-day intervals and are given in Table IV. 


TABLE IV—TABULATION OF (M4 - Mpg) FROM RITCHIE’S CURVE 


Daysof Oct. Nov. Dec. Feb. Mar. 
month 

I-10 
11-20 
21-31 


Jan. 
degrees Celsius 
+ 2*7 
+16 
+1°5 


Apr. 
+o°8 


+10 


+12 


+11 
+09 
+0°7 


+1°3 
+14 
+15 


+14 
+13 
+ 1°2 


+0°5 
+02 
0-0 


+ 2°97 
+1°7 
+17 


This table was then used to obtain values of M,4-— Me for two winters at 
Watnall and three at Wyton and the errors calculated. They are summarized 
in Table V. As Student’s é-test gives strong grounds for believing that the 
data for the two places are samples of the same population, the averages 
are also given at the foot of the table. These averages are very similar to the 
only comparable figures available from the use of Parrey’s method, namely 
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TABLE V—ERRORS IN (M,4—- Mp) BY RITCHIE’S METHOD 


Place Number of Total number Mean error S.D. 
winters of occasions 

degrees Celsius 

Wyton 3 560 — 0°09 1-08 

Watnall 2 279 +008 1:09 

Mean 0-00 1-09 


those for Wyton (see page 352). When the two sets of figures were compared 
by means of the f-test, a value of ¢ = 0-7 was obtained. As the critical value 
of ¢ at the 5 per cent level is 1-96, it may be accepted that both sets of figures 
refer to samples of the same population. This was to be expected as Parrey’s 
equation and Ritchie’s curve are simply different ways of expressing the 
dependence of M4 —- Mr upon the length of the cooling period. 


Combined with the data for errors in forecasting My, given by Steele, 
Stroud and Virgo (mean error — 1:10 degC, co, 2-79 degC) the figures for the 
whole process for Wyton are : 


mean error -— 1-09 degC 
standard deviation 3-00 degC. 


For Watnall for 1969-70 the mean error in forecasting My, is — 0-21 degC 
and the standard deviation is 1-99 degC, irrespective of whether occasions 
when D,, was greater than 10°C are included or not. The mean error of 
Ma -— Mp was zero with o 1-09 degC so that for the. whole process for Watnall 
the figures are : 


mean error -— 0-21 degC 
standard deviation 2-27 degC. 


Discussion. To clarify the discussion the essential figures have been 
gathered together in Tables VI and VII. For Wyton the direct regression 


TABLE VI—COMPARISON OF ERRORS IN FORECASTING Mr 
Method Source of observations Mean error c 
degC 


Direct regression with Wyton and Watnall 
cloud correction and combined +0-38 2°42 
rejecting D,,> 10°C 
Parrey’s method Wyton — 1°14 2°97 
Ritchie’s method Wyton — 1-09 3°00 
Ritchie’s method Watnall — O21 2°27 


TABLE VII—COMPARISON OF ERRORS IN FORECASTING (M4 - Mp) 


Method Source of observations Mean error S.D. 
degC 
Parrey’s Wyton ~ 0°04 1°02 
Ritchie’s Wyton and Watnall 0-00 1-09 
combined 


method gives better results, possibly because the large errors in forecasting 
Ma, at Wyton are masked to some extent in the direct regression equation. 
At Watnall, however, where the errors in forecasting M4 are much lower, 
Parrey’s and Ritchie’s methods give better results. For forecasting M4—- Mr 
it is immaterial whether Parrey’s or Ritchie’s method is used. 
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The agreement between the results for the two roads over 60 miles apart is 
remarkable and while it would be unwise to conclude that the tabulations 
used in these two roads are of widespread application, it may be assumed 
that these results show the order of accuracy attainable, and it is hoped that 
other investigators will be encouraged to try similar experiments on roads 
elsewhere. 
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551-509-323 
MAXIMUM TEMPERATURE ON CLEAR DAYS 


By G. A. INGLIS 


Summary. An account is given of some further developments of Gold’s technique for fore- 
casting maximum day temperature. The isothermal method is a variation of Johnston’s method 
in which the geometrical work is reduced; in the oblique method a set of el oblique straight 
lines is substituted for the isothermals; the saturated-adiabatic method makes use of a ‘mean’ 
saturated adiabatic instead of a ‘mean’ isothermal. In the thickness method the maximum 
temperature in summer is obtained from a table of 1000-850-mb thickness values. 


The isothermal method. Gold’s' method of estimating maximum day 
temperatures was modified by Johnston? and this in turn was slightly amended 
following a suggestion by Jefferson.* Johnston made use of assessments 
(denoted here by p,) made by Gold of the depth of the layer which is changed 
from an isothermal to a dry adiabatic by solar heating on clear days (Table I). 
The term ‘depth’ is used instead of ‘thickness’ which is used in a different 
sense later on. The curve does not really become a dry adiabatic, the lapse 
rate next to the surface being much steeper in the typical case, but Gold’s 
figures give good results in practice. As in all the tephigram methods the 
ascent used is that for 00 Gr. 


TABLE I—DEPTH OF LAYER (f,) WHICH IS CHANGED BY HEATING FROM AN 
ISOTHERMAL TO A DRY ADIABATIC, AND ASSOCIATED RISE IN TEMPERATURE (1,) 


Jan. Feb. Mar. Apr. May June July Aug. Sept. Oct. Nov. Dec. 
p,(mb) 60 80 95 110 120 125 120 4110 100 85 60 50 
7, (degC) 5°53 % £35 wh UU eS 9 7 5 45 


In Johnston’s method (Figure 1) the isobars py and fo-p, are drawn, 
corresponding to the pressures at the surface and the top of the layer. The 
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isothermal II’ is drawn to make equal areas between it and the environ- 
ment curve, IE the boundary of the upper area being a dry adiabatic. The 
point F where the dry adiabatic meets the surface isobar gives the forecast 
maximum temperature. 


~ 
Po eA 
Pi 








FIGURE I —JOHNSTON’S METHOD 


Gold also made assessments of the rise in temperature (r,) corresponding 
to the change from an isothermal to a dry adiabatic; in Table I these have 
been converted to Celsius to the nearest 0-5 degC. The values of p, and r,, 
can easily be verified by drawing on the tephigram a set of triangles FII’ 
for each month. 


If we know the value of r, it is not necessary to draw IF. If instead we 
note the temperature at I’ we can add to this the appropriate rise from 
Table I. This gives an alternative method in which part of the geometry is 
replaced by arithmetic, and it is considered to be somewhat easier and more 
accurate. 


Procedure. Draw the isobars py and p, — p, (Figure 2). Draw the isothermal 
II’ to make equal areas between it and the environment curve, IE the 
boundary of the upper area being a dry adiabatic. To the temperature 
at I’ add the rise (r,) from Table I to get the forecast maximum temperature. 
As in Johnston’s method it is helpful to use a transparent scale with two lines 
at right angles. 

Johnston says that the midnight environment curve is to be modified so 
as to approximate to conditions at dawn. But since modification of the curve 
is difficult and the change in area is usually insignificant it may be better to 
leave the curve unaltered. 


The oblique method. Johnston’s method has two advantages over Gold’s; 
it is independent of the scale of the tephigram and the areas to be estimated 
are relatively small because the midnight environment curve is often nearly 
isothermal. But this is not always so; in fact on the average the curve is 
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FIGURE 2—-THE ISOTHERMAL METHOD 


closer to the saturated adiabatic and often nearly coincides with it. It would 
be useful to have a technique based on estimating a ‘mean’ saturated adiabatic 
instead of a ‘mean’ isothermal. It was found that Johnston’s technique for 
the isothermals can be adapted to any other set of lines; but the saturated 
adiabatics are not suitable for use with a scale because they vary in slope. 
(The curvature is negligible.) Therefore another method — the oblique 
method — was devised which makes use of a set of parallel straight lines 
inclined at an angle to the isothermals. 

It is necessary to make assessments of the depth of the layer which is changed 
by heating from an oblique to a dry adiabatic and the corresponding rise in 
temperature. This could be done graphically as in Figure 3. For any month 
let I’F be the mean surface isobar (say 1020 mb), and let F correspond to the 


O | 





F 





FIGURE 3—GRAPHICAL METHOD OF ASSESSING VALUES OF fq, 1g, Ps, 13. AND 
DIAGRAM FOR OBTAINING FORMULAE FOR f, AND f 


mean daily maximum temperature for, say, southern England. Subtract r, 
from the maximum temperature to get the point I’ and let the isothermal 
and the dry adiabatic meet at I. Draw the oblique OO’ at the required 
angle of inclination « so that ORI and O’RI’ are equal in area. Then FO’ 
will give the rise in temperature r,, and the difference in pressure between 
O and the surface will give the depth of the layer py. 
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It is easier however and just as accurate in practice to calculate the values. 
If the surface isobar is assumed to be a straight line inclined at an angle of 
45° to the isothermals then by geometry 

Ts O’F I 


ry ——— /(1 + tan a) 
But r,p,; = rg, since area I1’F = area OOF, therefore 


Po 
Py 
These formulae were used to calculate, for an angle of inclination « of 45°, 


the values in Table II, and the values were then used to make the transparent 
scale illustrated in Figure 4. (not exact). 





=~yY(1+tana) . 


TABLE II—DEPTH OF LAYER (f,) WHICH IS CHANGED BY HEATING FROM AN 
OBLIQUE AT 45° TO A DRY ADIABATIC, AND ASSOCIATED RISE IN TEMPERATURE (r2) 
Jan. Feb. Mar. Apr. May June July Aug. Sept. Oct. Nov. Dec. 


f,(mb) 85 115 35 155 170 175 #+%\1I170 155 > 120 85 70 
ry(degC) 3545 SF PR MCHiL PE ops MP gy oogse veg 





FIGURE 4—-SCALE FOR USE IN OBLIQUE METHOD 


Procedure. The scale is used as follows (Figure 5). Lay the thin line 
corresponding to the month along the surface isobar. Adjust the position 
of the heavy oblique line OO’ so that equal areas are enclosed between 
it and the environment curve, the upper and lower boundaries being the dry 
adiabatic and the surface isobar. Note the temperature at O’ and add r, 
to get the forecast maximum temperature, which can also be read at F. 


The dimensions of the scale depend on those of the tephigram. Figure 4 
is }? the scale suitable for the tephigram currently in use, i.e. Metform 2810. 
Alternatively a scale consisting simply of two lines OO’ and OF inclined 
at an angle of 45° can be used. The procedure would then be to subtract 
p, from the surface pressure and adjust the position of OO’ so that the point 
O lies on the isobar p,-~, and OO’ intercepts equal areas between it and 
the environment curve as above. The temperature at O’ on the surface 
isobar is then read and the value of r, added. 
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FIGURE 5—-METHOD OF USING SCALE 


The saturated-adiabatic method. This method is for use without a 
scale. The necessary table of assessments of values of depth and rise in 
temperature could be made graphically by a construction similar to that in 
Figure 3, OO’ being in this case replaced by a saturated adiabatic SS’. But 
since the curvature of the saturated adiabatics is negligible they can be treated 
as straight lines, so that the method becomes a modification of the oblique 
method, in which the slope of the oblique varies from month to month. The 


values of the angle of slope in Table III were based on mean monthly values 
of the 1000-850-mb thickness for south-east England but a similar result 
can be got by drawing on the tephigram a set of saturated adiabatics through 
the points corresponding to the mean daily surface temperatures. From 
these values Table IV was calculated using the formulae for the oblique 
method. 


TABLE III—MEAN INCLINATION TO THE ISOTHERMAL OF THE SATURATED ADIABATICS 
IN THE ‘HEATED LAYER’ 


Jan. Feb. Mar. Apr. May June July Aug. Sept. Oct. Nov. Dec. 
angular degrees 


58 58 57 55 5! 47 45 45 49 5! 55 57 


TABLE IV—DEPTH OF LAYER (f,) WHICH IS CHANGED BY HEATING FROM A 
SATURATED ADIABATIC TO A DRY ADIABATIC, AND ASSOCIATED RISE IN TEMPERA- 
TURE (r3) 

Jan. Feb. Mar. Apr. May June July Aug. Sept. Oct. Nov. Dec. 
f3(mb) 100 190 «©6150 «€©6170)«©6180) = 180) 170s 15514525 95 80 
rs(degC) 3 4 5 6 7 ve. oe 7 6 5 3 3 


Procedure. Draw the surface isobar p, and the upper isobar fy — p, (Figure 6). 
Draw a saturated adiabatic SS’ making equal intercepts with the environment 
curve, SE the boundary of the upper area being a dry adiabatic. To the 
temperature at S’ add r, to get the forecast maximum temperature. The 
latter can also be read at the intersection of the dry adiabatic and surface 
isobar. In drawing the saturated adiabatic the curvature can be disregarded. 





Meteorological Magazine, 99, 1970 


FIGURE 6—THE SATURATED ADIABATIC METHOD 


The 1000-850-mb thickness method. In the methods described above 
we see that the maximum temperature can be regarded as consisting of two 
parts — the ‘mean’ temperature in a layer of the atmosphere next to the 
surface, which is a characteristic of the air mass, plus a constant which depends 
on the month. But the mean temperature of the layer between two pressure 
levels can also be measured by its thickness in geopotential metres, and in 
the summer months it happens that the layer examined to find the Gold 
maximum is very nearly the same as the 1000-850-mb layer at 00 Gar. 


The areas equated in the isothermal method are not exactly the same as 
those equated in calculating the 1000-850-mb thickness; if they were the 
same then the thickness value, from a suitable table, would give the Gold 
maximum exactly. Actually two curves with different profiles but with the 
same thickness and same surface pressure have Gold maxima which differ by 
an amount which is usually less than 1 degC. If, in order to make a suitable 
table, we assume an average value of the lapse rate then the thickness value 
can be used to give an estimate of maximum temperature very close to the 
Gold maximum. 


The average value of the lapse rate was assumed to be three-quarters of the 
saturated adiabatic and from the Radio Sounding Diagram (Metform 2813A) 
Table V was prepared, an allowance being made for an assumed mean relative 


TABLE V—VALUES OF 1000-850-mb THICKNESS AND ASSOCIATED SURFACE 
TEMPERATURE (AT 1020 mb) ASSUMING THREE-QUARTERS OF SATURATED 
ADIABATIC LAPSE RATE AND 75 PER CENT RELATIVE HUMIDITY 

Phickness 

gpm 1280 1290 1300 1310 1320 1330 1340 1350 1360 1370 1380 1390 
[emperature 

(°C) o-5 «614 3°3 52 71 gO 109 12°77 1475 163 18:2 20°0 


humidity of 75 per cent. (The tephigram is not suitable for calculating values 
of the 1000-850-mb thickness.) Table VI was then prepared by proportional 
parts from Tables I and IV, and Tables V and VI were combined in Table VII. 
To avoid the need for interpolation the maximum temperatures corresponding 
to each geopotential metre have been worked out for August in Table VIII, 
which can also be used for the other months if a small correction is applied. 
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TABLE VI—RISE IN TEMPERATURE ASSOCIATED WITH CHANGE BY HEATING FROM 
A CURVE WITH THREE-QUARTERS OF THE SATURATED ADIABATIC LAPSE RATE 
TO A DRY ADIABATIC 

April May June July August September 

Rise (degC) 6-9 8-0 8+5 8-3 78 6-8 

TABLE VII—VALUES OF GOLD MAXIMA FROM 1000-850-mb THICKNESS, ASSUMING 

THREE-QUARTERS OF SATURATED ADIABATIC LAPSE RATE AND 75 PER CENT 

RELATIVE HUMIDITY 

Thickness April ] June July August September 
gpm degrees Celsius 
1280 6-4 . 8-0 78 3 63 
1290 8-3 Q° 9°9 9°7 2 8-2 
1300 10°2 6 11°8 11°6 Lisl 10°1 
1310 12°1 “* 13°7 13°5 13°09 12:0 
1320 14:0 . 156 15°4 49 13°9 
1330 15'9 . 17°5 17°3 168 158 
1340 17°8 . 19°4 19°2 18-7 17°7 
1350 19°6 . 212 21-0 20°5 19°5 
1360 21°4 . 23°0 228 22°3 21°3 
1370 23°2 : 24°8 24°6 i 2371 
1380 251 , 26°7 26°5 26-0 25°0 
1390 269 28: 28:5 28-3 27:8 26:8 

TABLE Vvili— 1000-850-mb THICKNESS VALUES AND CORRESPONDING GOLD 

MAXIMA FOR AUGUST 

Thickness o 4 5 6 7 8 9 
gpm degrees Celsius 
1280 7°3 7 , . 8-1 8-3 8-4 8-6 8-8 
1290 g°2 9° , . 100 =: 10" 103 105 10°7 

1300 11] 11° or: ee a | 122 «6124 «= 128 

1310 130) 3=—s-« 13" 


8 * . . . 
20149 151 153155 157 159 100 1e2 164 
16° 


1330 17°0 . . 17°6 17°7 17°9 18-1 18-3 
1340 18-7 18-9 , . 194 196 198 200 20-1 
1350 20°5 20°7 ° . 202 214 8606216 )—atBeatg 
1360 22°3 225 ° 228. 230 232 234 ®©2«©2g6 8 23-7 
1370 24:1 9-243 , 24°7 + 24°9 a5" 252 254 25°06 
1380 26:0 =. 262 . 205 26:7 26:9 2771 27°30 —-27°4 
1390 278 280 282 2839 285 28:7 289 2g1 29:2 
April May June July September 
Correction (degC) -—0'9 +02 +0°7 +0°5 - 10 
The meaning of Table VIII is that, for example, in the month of August a 
curve with a thickness of 1335 gpm, a surface pressure of 1020 mb, a relative 
humidity of 75 per cent and a lapse rate three-quarters of the saturated 
adiabatic will have a Gold maximum of 17-7°C exactly. If with the same 
surface pressure and relative humidity the lapse rate is really isothermal, 
then Table VIII will overestimate the Gold maximum by an amount between 
o-7 and 1-0 degC. If the lapse rate is really a saturated adiabatic then the 
table will underestimate the Gold maximum by about 0-2 degC. 


Differences also arise from variations in surface pressure; 1020 mb was 
taken as the mean surface pressure on sunny days (i.e. slightly higher than 
the mean pressure for all days), and strictly speaking, a correction of about 
o-1 degC should be added or subtracted for every 2 mb above or below 
1020 mb; but this correction is usually small and was disregarded in the 
test. Other differences may arise from variations in humidity and from 
irregularities in the shape of the curve near 850 mb. 
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Test. A test of the method was carried out using data for the months 
April to September 1966-69. Days with 40 per cent or more of possible 
sunshine and no precipitation between og and 15 GMT were selected. The 
upper air ascent for 00 cmt at Aughton was used to forecast the maximum 
temperature by (i) one of the tephigram methods and (ii) the thickness method. 
(Aughton was used because plotted ascents from southern England were not 
available.) The forecast maxima were compared with the reported maxima 
at one or more of the stations Ringway (334), Shawbury (414), Birmingham 
(534) and Watnall (354). If more than one of these qualified as ‘sunny’ the 
mean was taken. 


TABLE IX-——-COMPARISON OF TEPHIGRAM AND 1000—500-mb THICKNESS METHODS 
AS PREDICTORS OF MAXIMUM TEMPERATURE USING DATA FOR AUGHTON 1966-69 


Number of Root-mean-square error 
observations Tephigram method Thickness method 
degC degC 


April 3: 167 1°56 
May 1°39 I° 

June 4 1°65 ee 
July 1°94 1-82 
August ‘ 1°55 1-68 
September 1°13 1-10 
Apr.—Sept. ‘ 1-62 1°59 


The results are shown in Table IX. It will be observed that there is no 
difference in accuracy between the two methods and hence in the summer 
months the thickness method could be used instead of the tephigram methods 
with equally good results. Although the method was fully tested only for the 
six summer months it is thought that it may be useful for some other months 
such as March and October. All that is necessary is to add these months 
to the list of corrections in Table VIII. 


The thickness method has several advantages. Firstly, there is no need to 
plot and examine the tephigram; secondly, a central office can obtain a 
forecast very quickly for a large area, e.g. north-west Europe, by plotting the 
thickness values on a chart; thirdly, the method can be used by inexperienced 
staff as it simply involves a subtraction followed by reference to a table. 


Regression technique. The direct way to establish a relationship between 
thickness and maximum temperature would be by means of a regression 
equation, using the same type of data as in the test above. It would have 
to be done separately for each month, except that the midsummer months 
could be grouped together. The result would be entirely independent of 
Gold’s figures and might be the best method if sufficiently large samples were 
taken. 


Comments. During the past two years the tephigram methods described 
above have been found to be easy and accurate. The curve is usually nearer 
to the saturated adiabatic than to the isothermal and is often so close to it 
that the areas involved are almost nil in which case the estimation can be 
done by eye alone. The reduction in geometrical work makes the methods 
more suitable for use with FAX copies of the tephigram than the older 
methods. The thickness method provides a quick and accurate estimate of 
the maximum temperature without reference to the tephigram. 
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551-509-324-2 551-557 
DEW-POINT TEMPERATURE AS A SNOW PREDICTOR 


By B. J. BOOTH 


Summary. Winter observations made at Upavon were analysed to see if there is any correlation 
between the dew-point temperature and the type of precipitation occurring at the time of 
the observation, which might allow the dew-point to be used as a snow predictor. A relationship 
was found and an empirical formula is suggested for use by the practising forecaster. 


Introduction. Many, if not all, of the snow predictors now in use resort 
to the use of upper air data. The various methods have been summarized 
by Boyden.? 

The use of upper air data creates three basic problems: 

(i) Which ascent or ascents are applicable to a particular area? 


(ii) The ascents are only available at 12-hourly intervals, and may 
substantially change in any 12-hour period. 


(iii) Most of the radiosonde stations are close to the sea, and in the lower 
levels may not be representative of a particular air mass. 


Boyden in his paper finds that the surface air temperature is an unreliable 
snow predictor. What is needed is an element which is readily available 
hour by hour, easily measured, and reasonably conservative. With this in 
mind an analysis was made of the full synoptic observations at Upavon, to 
see if it would be feasible to use the dew-point temperature as a snow predictor. 
Upavon is an exposed station some 580 ft (177 m) above mean sea level, 
and about 35 nautical miles (65 km) from the south coast. Full synoptic 
observations are made at 0g, 12, 15 and 21 GMT. Temperature observations 
are not available in the form of a continuous record. 


The data. The dry-bulb temperature, dew-point temperature and type 
of precipitation were extracted from the Upavon observations on all occasions 
when precipitation was occurring at the time of the observation, or had 
occurred during the previous hour, for the winter months of the years 1950 
to 1969 inclusive. No note was made of the intensity or duration of the 
precipitation. Winter in this case was defined as the months December, 
January and February. 


Results. An analysis was made of the dry-bulb temperature and 
associated precipitation so that a comparison could be made of the methods 
of forecasting snow using dry-bulb temperature and dew-point temperature. 
Scatter diagrams were then prepared, in which the type of precipitation was 
plotted against (i) the dew-point temperature and (ii) the dry-bulb temperature. 


It was noted that snow never fell with a dew-point >3°C and rain never 
occurred with a dew-point <- 3°C. 
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In the dew-point diagram the type of precipitation showed a marked 
change from rain to snow at about o°C. Boundaries were then marked on 
the diagram such that 2°C embraced values between 2-4°C and 1-6°C 
inclusive, 1°C embraced 1°5°C to 0°5°C, etc., ie. the approximations used 
when coding synoptic reports. The percentage numbers of occasions on 
which snow fell for each temperature were then calculated. Snow was 
classified as all types of freezing precipitation (excluding hail and sleet*). 
Occasions of sleet were also noted and are shown in the histograms, Figures 
1(a), 1(6), 2(a), 2(6). Freezing rain or drizzle was also included as snow, 
but this only occurred on g occasions in the 20 years. 


Total number of cases 
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Dew-point temperature (°C ) Dry-bulb temperature (°C) 
(a) With varying dew-point temperature. (6) With varying dry-bulb temperature. 
FIGURE I-—-HISTOGRAMS SHOWING THE FREQUENCY OF OCCURRENCE OF DIFFERENT 


TYPES OF NON-SHOWERY PRECIPITATION 
* Sleet : see footnote. 


The histogram in Figure 1(a) shows a marked increase in the frequency 
of snow when the dew-point changes from 1°C to o°C. The solid black area 
of each column shows the percentage number of occurrences of sleet. The 
total number of occasions examined was 438. 


An analysis of the dry-bulb temperatures and type of precipitation gave 
the histogram in Figure 1(6). The increase of frequency of snow occurrences 
is a much steadier process. Total number of cases examined was 406. 


Separate diagrams were produced for occasions when showers occurred. 
Although the size of the sample was small, the results, as shown by the 
histograms, Figures 2(a) and 2(b), show a similar pattern to the results 
previously noted. Here however the increased frequency of snow showers 
occurs with a dew-point of minus 1°C. Smith,’ in a recent paper, also noted 
that a dew-point of o°C was the critical value for the change from rain to 
snow. 

* In the United Kingdom the term sleet is used to denote precipitation of snow and rain (or 
drizzle) together or of snow melting as it falls, but it has no agreed international meaning. 
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Total number of cases 
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(a) With varying dew-point temperature. (6) With varying dry-bulb temperature. 


FIGURE 2——HISTOGRAMS SHOWING THE FREQUENCY OF OCCURRENCE OF DIFFERENT 
TYPES OF SHOWERY PRECIPITATION 


* Sleet : see footnote on page 364. 


As there seemed to be a marked increase in the frequency of non-showery 
snow with a dew-point of o°C a further analysis was made to see if a more 
accurate forecast of snow or rain/drizzle could be made, taking into account 
the dew-point depression. A scatter diagram was made with the dew-point 
as the abscissa and the air temperature as the ordinate. The type of precip- 
itation was plotted at the intersection of each value of dew-point and air 
temperature. The results are summarized in Tables I, II and ITI. 


TABLE I—SLEET* : FREQUENCY OF OCCURRENCES 


Dew-point Dew-point in degrees Celsius 
depression -2 -1 o I 
degC per cent 

3 16 

2 4 40 

I 4 25 15 

o 3 5 13 

* Sleet : see footnote on page 364. 


TABLE II-—-RAIN : FREQUENCY OF OCCURRENCES 
Dew-point Dew-point in degrees Celsius 
depression -1 °o I 2 

degC per cent 
3 
2 20 
I 29 
o 18 
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TABLE III—SNOW : FREQUENCY OF OCCURRENCES 
Dew-point Dew-point in degrees Celsius 
depression -2 -1 o 1 
degC per cent 
3 100(18) 51(6) 
2 88(8) 85(26) 40(5) 
i 100(24) 88(47) 46(28) 8(26) 
o = 97(38) 77(56) 15(54) 3(30) 
Numbers in brackets represent total number of occasions examined for each 
dew-point and dew-point depression. 


It can be seen from these tables that (i) snow is unlikely with a dew-point 

> 1°C; (ii) there is a good chance (> 50 per cent) of snow with negative 
dew-points even though the air temperature may be positive; (iii) there is a 
moderate chance (30-50 per cent) of snow with a dew-point of o°C even 
with a dew-point depression of 1 degC. 


The process was repeated for showery situations. Unfortunately the sample 
was not large enough to produce any conclusive results. 


Examination of Table III suggests that there is a relationship between 
dew-points < o°C (DP) and dew-point depressions (DPD) which gives the 
boundary between snow and precipitation other than snow. This is : 


DPD = [(-2xDP) + 1] degC, 
e. if DPD < [(-2 x DP) + 1] degC, then there is a good (> 50 per cent) 


chance of snow. 


Table IV gives the values of dew-point and dew-point depression required 
for > 50 per cent chance of snow using the above relationship. These values 
are easily obtainable from hourly synoptic reports. 


TABLE IV-——-VALUES OF DEW-POINT AND DEW-POINT DEPRESSION FOR PROBABILITY 
OF SNOW > 50 PER CENT 
Dew-point in degrees Celsius 0 -1 -2 -3 -4 
Dew-point depression in degrees 
Celsius equal to or less than 1 3 5 7 9 


Result of check on independent data. Using Table IV and adding 
the rider ‘no snow with dew-points 0°C’, the method was tested on independent 
data for March 1950-69. Although the frequency of snow and sleet is declining 
in March, it is still, at Upavon, as high as in December. The results are 
summarized in Table V(a). 


TABLE V(a)—-NON-SHOWERY PRECIPITATION 


Forecasts of rain Correct forecasts Forecasts of snow _ Correct forecasts 
49 37(76 we cent) 46 41(89 per cent) 


5 of the incorrect forecasts were to Bis 1 of the incorrect forecasts was due to sleet. 
TABLE V(b)—SHOWERY PRECIPITATION 
Forecasts of rain Correct forecasts Forecasts of snow Correct forecasts 
18 14(78 per cent) 25 23(92 per cent) 
Table V(6) gives the results for showery precipitation. Here, however, 
the criterion used was simply a forecast of snow with a dew-point < - 1°C. 


This cruder criterion was used mainly because of the lack of conclusive data 
as previously mentioned. 
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Discussion. As the analysis of dew-points and dew-point depressions 
was made without reference to the intensity or duration of precipitation, 
no conclusion could be made about the rate of increase or decrease of the 
dew-point during precipitation. One would expect a change of dew-point 
with the approach of a different air mass, and it is possible to follow the 
change of dew-points and dew-point depressions on a normal working chart. 
It is then possible to decide, particularly in the case of an approaching warm 
front in winter, when precipitation falling as snow will turn to rain. 


Lumb’ states that snow can extend down to the 1°5°C wet-bulb level; 
unfortunately surface wet-bulb temperatures are not plotted on synoptic 
charts. If the figures in Table IV are plotted on a tephigram and the wet- 
bulb temperatures calculated, it will be seen that they lie between 0-5°C 
and 1°5°C. Thus by using the criteria in Table IV we now have a method 
of estimating, hour by hour, from the synoptic charts the change in wet-bulb 
temperature at the surface. Lumb‘ shows how prolonged moderate or heavy 
precipitation can also cause a change of dew-point. 


The figures in Table IV were calculated for occasions when precipitation 
was occurring. It is felt, however, that in situations where precipitation is 
expected to occur, if the conditions of Table IV are fulfilled, then it would 
be reasonable to forecast the precipitation as snow. 


Table VI gives the 2 x 2 contingency table comparing forecasts of rain 
or snow with the actual occurrences of rain and snow for the test forecasts 
on the March data. The chi-square value (with Yates’s correction) is also 
shown. 


TABLE VI—NUMBER OF FORECASTS OF SNOW OR RAIN AT UPAVON COMPARED 
WITH THE TYPE OF PRECIPITATION THAT ACTUALLY OCCURRED 


Rain forecast Snow forecast Totals 
Rain occurred 37 5 42 
Snow occurred 12 4! 53 


Totals 49 46 95 
The chi-square value is 37-6, which is significant at the o-1 per cent level. 


The results of the test, Table V, are encouraging; however it would be 
interesting to see if similar results would be obtained from a low-level inland 
site, and a coastal site. 


The object of this paper was to see if there was a correlation between the 
type of precipitation and dew-point temperature; this would seem to be 
the case with the dew-point of o°C being the critical value for non-showery 
precipitation and — 1°C for showery precipitation. 
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551.558-21 
MOUNTAIN LEE WAVES IN THE VALE OF YORK 
By G. F. RUDDOCK 


Following publication of Casswell’s simplified method of forecasting lee- 
wave conditions,’ a number of aircraft reports were received at Linton-on- 
Ouse during frequent westerly spells in the period June 1966 to December 
1967. These westerly spells frequently produce wave conditions over the 
Vale of York, because of flow over the Pennines. 


Examination of the reports, together with ground observations of wave- 
type clouds, suggested that, broadly speaking, Casswell’s method worked 
fairly well in regard to maximum vertical velocities. However, the investiga- 
tion drew attention to two important variations from Casswell’s findings : 


(i) Casswell quoted 20 kt* as the likely lower limit of wind speed at 
the mountain crest for wave formation. In this area the incidence 
of wave clouds, or of waves reported by aircraft, is quite common 
with wind speeds down to about 15 kt. This is in agreement with the 
statement in WMO Technical Note No. 34, page 119, that the minimum 
speed for wave formation is about 7 m/s for small mountains. 


Casswell’s method involves smoothing the environment temperature 
curve between 1000 and 700 mb, and between 700 and 300 mb. 
This has the disadvantage of smoothing through inversions, which 
may be major discontinuities. Theory strongly points to the height 
of maximum vertical velocity being at the inversion level (see for 
example Corby,* para. 19). The examination of aircraft reports leads 
to the same conclusion. It follows that the height of the maximum 
vertical velocity, as deduced by Casswell’s method, is open to doubt 
if it differs widely from the base of the stable layer. 


A comparison was made for 57 occasions (all suitable for mountain waves) 
between the height of maximum vertical velocity as deduced from Casswell 
and the height of the base of the isothermal or inversion layer, as taken from 
the representative tephigram. 


This gave the following results : 


Number of occasions with base of stable layer below 10 000 ft* 49 
Number of occasions when Casswell method gave height of maximum 

vertical velocity within 1000 ft of the base of the stable layer 14 (28-6 per cent) 
Number of occasions when Casswell method gave height of maximum 

vertical velocity more than 1000 ft above the base of the stable layer 32 (65:3 per cent) 
Number of occasions when Casswell method gave height of maximum 

vertical velocity more than 1000 ft below the base of the stable layer 3 (6-1 per cent) 

If the maximum vertical velocity occurs at the base of the stable layer 
then the Casswell method will give too high a level on a substantial number 
of occasions presumably because of the smoothing used in the procedure. 
Thus the height of inversion, or of base of isothermal, should be used to fore- 
cast the height of maximum vertical velocity, rather than the theoretical 
value. In practice, the nearer the stable layer to the mountain top, the more 
important this becomes. Light aircraft with maximum power lift of 1000 to 





* ikt & 0-5 m/s; 1000 ft & 300 m; 1000 ft/min & 5 m/s. 
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1500 ft/min have insufficient power to overcome downdraughts of 1200 to 
1500 ft/min and these are not uncommon in waves due to the Pennines. 
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REVIEWS 


Atmospheric tides, by Sydney Chapman and Richard S. Lindzen. 240 mm x 
165 mm, pp. ix + 200, illus., D. Reidel Publishing Company, 419-421 
Singel, PO Box 17, Dordrecht, Holland, 1970. Price : Dfi. 38. 


This monograph must be one of the last of the great contributions made 
by Professor Sydney Chapman to our knowledge of the upper atmosphere. 
It was with deep regret that the members of a Symposium of the American 
Meteorological Society on “The Dynamics of the Mesosphere and Lower 
Thermosphere’ at Boulder, Colorado (Bull. Amer. met. Soc., Lancaster, Pa, 51, 1970, 
Pp. 377-387), learned that because of illness he would not be able to be 
Chairman of one of the sessions and later that he had died. He was 82. The 
first of his papers quoted in the references was dated 1913 and the last in 
1969 and there were several in each decade between. As all interested in 
atmospheric physics will know, his contributions were not restricted to the 
tides. Other subjects such as the ionosphere, the earth’s magnetism, ozone 
photochemistry, solar radiation, the theory of non-uniform gases, etc., etc., 
have equally benefited from his great insight. His co-author Professor Lindzen 
of the University of Chicago some 50 years his junior has already established 
himself as a world authority on tidal and other wave motions in the upper 
atmosphere. Their joint review of the present state of knowledge of this 
complicated subject is authoritative and must become a standard reference 
book for many years. 


Attempts to explain theoretically the magnitudes and phases of the observed 
tides in the ocean and in the atmosphere have been made by outstanding 
physicists and mathematicians for almost the last three centuries and a short 
discussion of the contributions due to such scientists as Newton, Laplace, 
Kelvin, Rayleigh, Lamb, Taylor and many others are included in the 
introductory first chapter which also gives a clear account of the fundamental 
concepts and outstanding problems. It is known that the gravitational pull 
of the moon is greater than that of the sun and hence in the ocean the lunar 
tides are greater than the solar. In the atmosphere, however, the lunar tide 
measured by barometers is comparatively small and hence it appears that 
the thermal drive due to the sun must play a dominant role. Observations 
show also that the semidiurnal barometric solar tide is unexpectedly great 
compared with the diurnal component although the thermal forcing for the 
latter is clearly the greater. Earlier explanations in terms of a natural resonance 
of the atmosphere have now been rejected in favour of increased efficiency 
of coupling between the semidiurnal forcing terms and their induced motions. 
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The second and third chapters (labelled 2S and 2L) deal respectively with 
the solar and the lunar daily atmospheric oscillations as revealed by meteor- 
ological data. Near the surface the details are not easy to extract. The 
necessary harmonic analyses require careful statistical techniques and the 
detailed results are described at some length. At higher levels the tidal 
components become successively larger until above 100 km they dominate 
the general circulation of the atmosphere. The rapidity of the wind changes 
and their magnitudes and phases at different times and locations raise severe 
practical problems as regards methods of measurement and _ theoretical 
questions regarding wave propagation, forcing functions and dissipation 
mechanisms. 


The final chapter (3) sets out in detail the present position on the quantita- 
tive theory of tides. The basis as in all dynamical studies of atmospheric 
motions rests on the equations of motion and the continuity and thermo- 
dynamic equations and the analytical solutions of their linearized forms with 
the tidal constraints are given in full. The subject, however, is complex and 
the Hough functions which arise in the solutions have in practice to be 
evaluated numerically using electronic computers. Recent computations 
are presented in some detail, compared with observations and outstanding 
problems discussed. 


The division of labour between the two authors was broadly that Professor 
Chapman wrote most of Chapters 1, 2S and 2L except the sections dealing 
specifically with recent upper atmosphere observations and modern theoretical 
calculations while Chapter 3 was written by Professor Lindzen. Even if it 
had not been stated in the Preface this could have been deduced directly. 
In his section Chapman sets out his procedures in great detail while Lindzen 
prefers to assume that his readers do not need quite so much guidance. This 
may be a disadvantage as many students may be deterred by the difficulty 
of the theoretical aspects in some sections, leaving these to the specialists. 
In particular it might have been preferable to include more detailed explana- 
tions in Chapter 3 possibly with the inclusion of examples of the procedures 
involved. 


The subject is still moving rapidly due, to a large extent, to further work 
by Lindzen himself and no doubt in a few years a second edition will be called 
for. The sections covered by Chapman, however, will need little revision 
and until that time this edition will remain the major basic text and statement 
of all that was known about this subject in early 1970. The book is attractively 
set out and clearly illustrated but a few minor errors were noted presumably 
due to insufficient proof reading. 


R. J. MURGATROYD, O.B.E. 


Statistische Auswertung geophysikalischer und meteorologischer Daten, by J. Tauben- 
heim. 230 mm X 165 mm, pp. 386, illus. Geest & Portig K.-G., 
Sternwartenstrasse 8, Leipzig, 1969, Price : DM 61. 


The application of statistical methods to geophysical data has been 
surprisingly unsystematic, considering the number, variety and quality of the 
data and there is a need for an authoritative and comprehensive account 
of the tools of the statistician’s trade which is sympathetic to the geophysical 
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medium. Professor Taubenheim’s book goes a long way towards satisfying 
this need by providing a well-balanced and well-illustrated account, first of 
the standard statistical concepts and methods and then of subjects such as 
time series analysis, power spectrum analysis and filtering which are relatively 
more important in geophysical work than they are elsewhere. The illustrations 
are nearly all drawn from meteorology or other geophysical sciences and 
the depth of the theoretical treatment is enough to serve the needs of most 
professional meteorologists without being violently excessive. The treatment 
does not generally extend to the statistical aspects of methods, such as principal 
component analysis or stepwise multiple regression, which can only be exploited 
by computer, nor to subjects such as the statistical control of errors in data, 
which have become increasingly important in the era of computerized 
research. 


Some statistical tests such as Kolnogoroff’s and Wilcoxon’s are included 
which are not usually found in comparable accounts, and the general 
impression is of clarity and thoroughness, and of the collection in one place 
of material which is hard to find elsewhere. The bibliography is good and 
the index excellent. Altogether a book to be recommended, and this reviewer 
at any rate hopes that it will appear before long in English translation. 

J. M. CRADDOCK 


_ Meteorological data catalogue. International Indian Ocean Expedition Meteorological 
Monographs, Number 3, by James R. Nicholson. 210 mm xX 280 mm, 
pp. 59, tllus., East-West Center Press, Honolulu, Hawaii, 1969. Price: $7.50. 


This monograph is the third in a series covering the meteorological work 
of the International Indian Ocean Expedition (IIOE) between 1959 and 
1965. The author, an Environmental Science Services Administration (ESSA) 
Weather Bureau research meteorologist assigned to the IIOE, spent the 
years 1963 and 1964 at the International Meteorological Centre, Bombay, 
followed by three years at the University of Hawaii where he supervised data 
handling for this series of monographs, and so can write on the subject with 
some authority. As implied by its title the monograph is primarily intended 
to serve as a reference catalogue of data collected in the course of the 
meteorological and oceanographic research programmes of this expedition, 
and the reader will soon appreciate how well this purpose has been achieved. 
In the field of meteorological telecommunications there is material concerned 
with organization and operational experience, which could have a bearing 
upon planning for similar large-scale projects in the future within the World 
Weather Watch. Problems met with in the recording, collection and process- 
ing of data are next briefly described, followed by an account of methods 
of archiving which includes tables giving full details of sorties and observations 
made by research aircraft over the Indian Ocean, together with a description 
of instrumentation available on the aircraft and time-lapse films taken on 
these flights. Another table, summarizing all computer programme: (title, 
language, purpose, etc.) already written in connection with processing 
numerous categories of data collected during the IIOE, could be a valuable 
aid to research workers in these fields. 


The monograph concludes with a brief discussion and, arising from 
experience gained by the author and his colleagues during the IIOE, some 
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suggestions for future procedures for data collection and handling during 
World Weather Watch. While some valid points are made relating to quality 
control of data and modification of codes, it seems doubtful whether difficulties 
involved in revising codes at the international level have been fully appreciated 
by the writer of this book. 


The monograph is lavishly produced and contains a useful set of up-to-date 
references. Not surprisingly it is rather expensive in relation to its slim 
dimensions. 

R. F. M. HAY 


Numerical analysis — the mathematics of computing, Volume 1, by W. A. Watson, 
T. Philipson and P. J. Oates. 227 mmx149 mm, pp. xi+224, tlus., 
Edward Arnold (Publishers) Ltd, Woodlands Park Avenue, Woodlands 
Park, Maidenhead, Berks., 1969. Price: 25s. (limp edition). 


Those of us who have learned something of numerical analysis as and when 
it became necessary will welcome the publication of this book. Although 
designed primarily as a textbook for the special G.C.E. A-level course in 
numerical analysis, this book will be useful also in other sixth-form courses 
and in courses of further education. It provides a basic introduction to 
numerical methods of solving mathematical problems, such as must surely 
now be regarded as an essential feature of scientific education in schools. 
An increasing use of books such as this in science sixth forms will do much, 
I think, to eliminate confusion and apprehension from attitudes to the use 
of computers in scientific work. 


The first three chapters of the book prepare the ground for the remaining 
six by introducing the reader to hand calculating machines, flow charts, 
and the sketching of simple functions respectively. I found that the emphasis 
on hand calculating machines in the opening chapter gave the book a rather 
dull beginning. The rather tedious details of register setting and crank 
rotating might perhaps have been better placed in an appendix to the book. 
By contrast, I felt that the second and third chapters provided clear and 
attractive introductions to their respective topics. 


The rest of the book deals with a series of basic numerical techniques which 
include iterative methods for the solution of equations, differences, linear 
interpolation, and numerical integration up to Simpson’s rule. Chapter 
six, dealing with the solution of linear simultaneous equations, is especially 
useful, I think, introducing as it does such topics as the method of relaxation 
in a painless way. 


Throughout the book clear indication is given of the limitations of the 
techniques described and of the rounding and other errors involved in various 
operations. Profitable use is made of worked examples to explain the 
numerical methods, and a good selection of examples for the student (with 
answers) is included. Taken as a whole, there is no doubt that this is a useful 
introductory textbook on numerical analysis, and I look forward to the 
publication of Volume 2 in due course. 

A. J. GADD 
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Numerical analysis — the mathematics of computing, Volume 2, by W. A. Watson, 
T. Philipson and P. J. Oates. 227 mm x149 mm, pp. x+166, illus., 
Edward Arnold (Publishers) Ltd, Woodlands Park Avenue, Woodlands 
Park, Maidenhead, Berks., 1969. Price: 28s. (limp edition). 


The first volume of this work gave a useful introduction to some basic 
concepts of numerical analysis and encouraged one to look forward to the 
publication of this second volume. In the event, however, I found the 
second volume rather disappointing. In general one would not quarrel very 
much with the presentation of individual techniques; the disappointment 
stems rather from the suspicion that the student who follows through these 
two volumes may be left with the impression that numerical analysis is a 
tedious and confusing subject. 


Admittedly the authors had no easy task of presentation in this second 
volume. They are concerned in part with extending and generalizing the 
material of Volume 1 and in doing so they inevitably become involved with 
more complex techniques which lack both mathematical elegance and obvious 
usefulness. Even so, I feel that a more attractive presentation could have 
been achieved if an overall direction of development of the subject were 
more apparent. Volume 2 ends with a short chapter on the summation of 
slowly convergent series which is a complete anticlimax. One is tempted 
to believe that the fact that the book has three authors may have contributed 
to a certain lack of coherence from one chapter to another. 


More than half of the pages of this second volume are devoted to material 
which in earlier books on numerical analysis has generally appeared in the 
chapter headed ‘interpolation’. I find it difficult to judge whether or not 
such an emphasis is justified, but certainly a rather sustained effort is 
required to remain interested in the various details presented. 


As the authors admit, the inclusion of an introductory chapter on the 
numerical solution of differential equations has the difficulty that the readers 
of the book in schools may not be in a position to cope with the methods of 
solution which are most widely used. Some simpler methods are therefore 
presented, in the hope that an enthusiasm for further study will be stimulated. 
The authors deal mainly with first-order equations, despite the fact that 
the simplest numerical process is that for a second-order equation with the 
first derivative absent. The description given of Fox and Goodwin’s method 
for first-order equations is, I think, unnecessarily daunting, and would have 
been better introduced by way of the trapezium rule for integration. 


One feature of the layout of the printed material deserves comment. 
The tables, formulae, and numerical examples are printed so boldly that 
they overshadow the explanatory text and tend to obscure the continuity 
of the development. I liked the occasional historical references and thought 
that they could have been usefully extended to give more perspective to the 
place of numerical analysis in the general pattern of mathematics and physics. 


I remain grateful that the attempt has been made to provide an attractive 
introduction to numerical analysis at sixth-form level, but I feel that there 
is scope for an improved presentation in the near future. 

A. J. GADD 
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Field guide to snow crystals, by E R. LaChapelle. 210 mm Xx 150 mm, pp. v + 
101, illus., University of Washington Press, Seattle and London, c/o 
American Universities Publishers Group, 27-29 Whitfield St, London, W1. 
1970. Price: 62s. 


This slim volume consists principally of a collection of 65 very beautiful 
photographs of snow crystals, each of which is accompanied by a short 
commentary describing the situation in which the crystals were found and the 
important physical mechanisms involved in their formation. The photographs 
have been carefully selected, not so much for their symmetry and beauty, but 
rather to illustrate clearly the various forms which snow crystals on the ground 
can assume. In this respect the new book differs significantly from the classic 
collections obtained by Bentley and Humphreys and by Nakaya, who were 
primarily concerned with the forms of snow crystals as they first fall from 
the sky. 


In the opening chapter, LaChapelle describes very clearly the processes 
which lead to changes in the structure of precipitated snow when it is lying 
on the ground. He also gives an account of a new system for classifying snow 
on the ground which has been proposed by R. A. Sommerfeld and himself. 
In this scheme three major processes leading to changes in the structure of 
snow are recognized—equitemperature (destructive) metamorphism, temp- 
erature-gradient (constructive) metamorphism and firnification (the process 
by which snow changes into glacier ice). 

The book contains a short section explaining how best to observe and 
photograph snow crystals with fairly simple equipment, so that the interested 
reader can experiment for himself. 


In short, this book fully lives up to the claim made in its title to be a field 
guide to snow crystals and it will make a handsome addition to the libraries 
of all who are interested in the fascinating forms of snow crystals. 

J. T. BARTLETT 


Tropical Indian Ocean clouds, International Indian Ocean Expedition Meteorological 
Monographs No. 4, by Andrew F. Bunker and Margaret Chaffee. 282 mm x 
222 mm, pp. 194, tllus., East-West Center Press, Honolulu, Hawaii, 1970. 
Price: $10.00. 


This monograph, the fourth in the series reporting data-gathering activities 
during the International Indian Ocean Expedition, presents schematic cloud 
patterns deduced from flights over the Indian Ocean during 1963 and 1964. 
Cloud data were derived from time-lapse film frames exposed from aircraft 
at intervals ranging from 2 to 10 seconds. A vast number of photographs 
were taken but only one frame in every hundred was used in this analysis. 
The photogrammetric technique employed to determine the distances of 
cloud from the aircraft, and the heights of bases and tops, is described. The 
volume is not, however, a collection of cloud photographs as the title might 
suggest. 


Two main sections deal with the south-west and north-east monsoons 
separately and each section is liberally illustrated with diagrams. There 
are only 20 pages of text compared with 165 pages of diagrams. There is the 
surprising total of 346 figures, nearly all of them line drawings; only 23 show 
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actual cloud photographs and of these only 4 were taken from aircraft, the 
remaining 19 being from satellites. All photographs are of rather poor quality. 


The schematic cloud distributions along flight paths are presented in 
cross-section or map form with the minimum of discussion, and are supported 
by daily wind maps showing streamline and isotach analyses for the surface 
and 700-mb levels. These maps, 210 of them, are very much smoothed and 
do not show the wind data upon which they were based. In some cases the 
analyses are significantly different from others for the same days published 
in Monograph No. 1 of this series. For example, Figure 125 should be compared 
with Figure 21 of Monograph No. 1, where basic data are reproduced on all 
the maps and small-scale distortions in the wind field are shown to exist. 
Such small-scale variations may be more closely related to cloud formations 
than the broad-scale flow itself, but the authors do not attempt any detailed 
interpretation of cloud patterns in relation to the wind field. 


No legend to the symbols used in depicting clouds appears with the 
diagrams. The reviewer found it a little distracting to refer back contin- 
ually to the description of the method of representation given on page 8, 
and even then some of the schematic forms remained difficult to interpret. 
Some radiation and temperature measurements are discussed briefly and 
presented as composite — not daily — data. In particular the temperature 
data are of little help in studies of the daily cloud structures. 


The volume is presented in the same attractive way as the others of the 
series. The reviewer would have preferred the daily wind-field maps to be 
adjacent to the appropriate cloud charts instead of being gathered together 
at the end of each section — but this is a minor criticism. The book is likely 
to have little general appeal because it is really intended for the research 
worker interested in the weather of the Indian Ocean, and to him it presents 
a wealth of information for study and digestion. 

J. FINDLATER 


HONOUR 
The following honour was announced in the Queen’s Birthday Honours 
List 1970 : 
B.E.M. 
Mr V. Efstathiou, Communications Superintendent, Main Meteorological 
Office, Episkopi. 
AWARDS 


WMO Research Prize 
We note with pleasure the decision of the Executive Committee of the World 
Meteorological Organization to award the WMO Prize for 1970 to the 
following young research workers: Mr F. B. A. Giwa (Nigeria), Mr M. 
Yamasaki (Japan), Mr P. E. Merilees (Canada) and Mr F. Bretherton 
(United Kingdom), for research work in various fields of meteorology. 


Fifteenth IMO Prize 
This prize, which is awarded each year as a token of remembrance of the 
International Meteorological Organization, which was the non-governmental 
body preceding WMO, was awarded posthumously to Professor R. Scherhag, 
Berlin, who died on 31 August 1970. 
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NOTES AND NEWS 
Meteorological Magazine : price increase 
As from January 1971, the price of an issue of the Meteorological Magazine 
will be 4s. [20 p] and the annual subscription will be £2 145s. [£2.70] including 
postage. 


CORRECTION 


Meteorological Magazine, September 1970, p. 270. In the second line of the 
summary and the fourth line of the main text of the article on a heated 
anemometer, for ‘Mount Olympus, Greece, read ‘Mount Olympus, Cyprus’. 
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